Prompt Multi-Gluon Production in High Energy Collisions 
from Singular Yang-Mills Solutions 
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We study non-perturbative parton-parton scattering in the Landau method using singular 0(3) 
symmetric solutions to the Euclidean Yang-Mills equations. These solutions combine instanton 
dynamics (tunneling) and overlap (transition) between incoming and vacuum fields. We derive 
a high-energy solution at small Euclidean times, and assess its susequent escape and decay into 
gluons in Minkowski space-time. We describe the spectrum of the outgoing gluons and show that it 
is related through a particular rescaling to the Yang-Mills sphaleron explosion studied earlier. We 
assess the number of incoming gluons in the same configuration, and argue that the observed scaling 
is in fact more general and describes the energy dependence of the spectra and multiplicities at all 
energies. Applications to hadron-hadron and nucleus- nucleus collisions are discussed elsewhere. 



I. INTRODUCTION 

Tunneling phenomena related with topology of Yang- 
Mills fields are described semiclassically by instantons 
K10. Some manifestations of these effects related with 
explicit breaking of U(l) and spontaneous breaking of 
SU(Nf) chiral symmetries are by now understood in sig- 
nificant detail due to strong ties to hadronic phenomenol- 
ogy and lattice studies, see e.g. [|| for a review. 

We know much less about the role of instantons in 
cross sections of various high-energy reactions. Such 
studies got considerable attention in the early 90's, when 
baryon-number violating processes through instantons in 
the standard model have been actively discussed 
These studies were extended to hard QCD processes with 
small-size instantons, and there are still attempts to see 
their effects on multi-jet production at HERA (for recent 
review see pi). 

Recently [|,|| it was suggested that nonperturbative 
configurations composed of instanton/antiinstanton play 
an important role in parton-parton scattering amplitudes 
at high-energy, and may account for most of the soft 
pomeron slope and intercept. The logarithmic rise of the 
inelastic cross-section was shown to follow from coher- 
ent multiple gluon production as described by the semi- 
classical field following from an interacting instanton- 
antiinstanton configuration. This mechanism was shown 
to be the same for pp and pp B: so no odderon ap- 
pears in the classical limit. At low energy transfer in the 
center of mass, the interaction is dipole, and accounts 
for the rise in the partial cross section from first prin- 
ciples. At intermediate and large energy transfer, the 
dipole approximation is not reliable as strong instanton- 
antiinstanton interactions set in together with unitarity 
constraints [||||. 

It was realized then that unlike perturbative processes, 
for which production of each subsequent gluon is as- 
sociated with a power of small coupling constant, in 
instanton-induced processes the cross section rises with 



a number of produced gluons reaching a maximum at 
some parametrically large value N « l/a(p). The phys- 
ical reason for this is that specific coherent clusters of 
gauge field are produced instead of independent gluons. 
In a recent paper [fl0[ the properties of minimal clusters 
of such kind were discussed: the clusters themselves were 
identified via constrained minimization of the Yang-Mills 
energy, with fixed size and Chern-Simon number, and 
their subsequent real time evolution have been studied 
both numerically and analytically. 

However, in real collisions the Yang-Mills energy or 
action of the final state is only one of the contributing 
factors. Another crucial factor is the overlap between 
the initial system of colliding gluons and the instanton, 
or whatever the tunneling path is. A particularly useful 
way jlffl to treat both semiclassical factors from first prin- 
ciples together, and enforce the unitarity constraints is to 
use a semiclassical approximation to the partial cross sec- 
tion based on an adaptation of the Landau formula for 
overlaping matrix elements in terms of singular field con- 
figurations. The occurence of a singularity is an essential 
feature of the Euclidean field configuration that interpo- 
lates between the vacuum at t = — 00 with zero energy, 
and the escape point at t = with finite energy. 

Diakonov and Petrov H used this approach (Landau 
method) to rederive the low-energy limit of the ampli- 
tude, and obtain a new result in the high-energy limit. 
Although they have not worked out the gauge configura- 
tion at small times explicitly, they were still able to as- 
sess the corresponding action, tunneling time and energy, 
to predict the behavior of the pertinent cross section at 
high-energy. A comparison of the low- and high-energy 
limits, has confirmed that the cross section has a maxi- 
mum near the sphaleron energy with a value that is close 
to the square-root of the cross section near zero-energy. 

The aim of the present paper is to analyze further the 
singular gauge configurations at the escape time t = 0. 
Following j(| we start with the high-energy limit and 
show how one can find the gauge configuration at the 
turning point which is the produced gluon cluster at 
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t = 0. It turns out that this configuration can be related 
to minimal YM sphalerons fic|| by a particular scaling 
law, containing a power of Q/M s . We show that the sub- 
sequent evolution (explosion) in real time of these clusters 
can also be derived, thereby generalizing recent work ]Tc| ] 
by our rescaling. The resulting spectrum of gluons and 
their multiplicity are readily obtained. We argue that the 
prescribed scaling law is valid not only at high energies, 
above the sphaleron mass, but in fact at all energies. 

In section 2 we introduce the standard notations for 
the spherically symmetric gauge configurations used. In 
section 3, we recall the key points in assessing the inelas- 
tic parton-parton scattering cross-section in the eikonal 
approximation. In section 4, we show that the singu- 
lar gauge configurations at the escape points follow from 
the sphaleron point for all energies through a pertinent 
rescaling, which is our main result. In sections 5, we 
assess the number of incoming and outgoing gluons per 
sphaleron for fixed center of mass energy in the semiclas- 
sical approximation. Our conclusions are in section 6. 
The Appendices contain a number of useful results in- 
cluding a stability analysis of the escaping configurations 
under perturbative light pair-decay. 



II. 0(3) SYMMETRIC YANG-MILLS 

We consider the QCD Yang-Mills theory wherein all 
dimensions are rescaled away by the sphaleron (anti- 
sphaleron) mass Ms, 



M s = 



1 

4a 



dx r 2 



96p 4 



3tt 



(x 2 + p 2 ) 4 Aap 



(1) 



unless specified otherwise. In the vacuum a ~ 0.3, p ~ 
1/3 fm with typically Ms ~ 3 GeV. In the scattering 
process, the sphaleron (antisphaleron) size p may change. 
We work mainly in temporal gauge. The YM potential 
V, kinetic energy K and Chern-Simons number N are 



V = 



K 



M s 
47r a 
M s 
47r a 



J dx±(F t 



a \2 
ij I 



N = 



^ J dxe ijk (A1 d 3 A a k + l -e abc At A) A%) 



(2) 



where we have ignored quarks for simplicity. 

In parton-parton scattering with large y/s, the incom- 
ing kinematics boils down to an eikonalized cross sec- 
tion as in (||) with a partial cross section a (Q) with 
Q w Ms <C \fs. In the center of mass it is reasonable to 
consider the gluonic configurations that maximizes a (Q) 
to possess higher symmetries than completely arbitrary 
fields. Here we take them to have spherical 0(3) symme- 
try of the sphaleron 



, t) = +e aij n, (1 - A(x, t))/x 
+ (5ai - n a rn) B(x, t)/x 
+n a riiC(x,t)/x 



(3) 



where rij = x*i/x is a unit 3- vector, and \x\ — x > a 
radial variable. Since we arc interested in singular gauge 
fields, we assume the three 2-dimensional independent 
functions A, B, C to be continuous and differentiable ev- 
erywhere in Euclidean space, except at x = where a sin- 
gularity will be located for fixed times ±T/2. In terms 
of (||) the Euclidean action 



S = dt(K + V) 



(4) 



reads 
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- C 2 + A' 2 + B' 2 



1 B 2 + 
C 2 {A 2 +B 2 ) 2C{A'B-AB') 



2x 2 



) (5) 



where the time variable has been rescaled through 
tMs — > t, in (||). The integration interval is to be spec- 
ified below in the presence of a time-singularity. Note 
that the energy 



= K V 



(6) 



where in Euclidean space —V plays the role of the po- 
tential. For self-dual configurations Q = 0. The Chern- 
Simons number is 



N = — / dx I A'B — B'A + -(A 2 + B 2 - V 
2tt J \ x 



(7) 



The dual to the 0(3) ansatz used here that maps on 
the antisphaleron follows similar reasoning and results as 
shown in Appendix A. The difference is a negative N. 



III. THE EVOLVING VIEWS ON 
NON-PERTURBATIVE HIGH ENERGY 
PROCESSES 

In this work we discuss semi-hard parton-parton pro- 
cesses at fixed ~ 1 GeV <C x^s, using nonpertur- 
bative QCD gauge configurations related to topological 
tunneling. The original approach [Q wich we followed 
in H was based on semiclassical instanton solution in 
the amplitude. To leading order in the small tunneling 
or diluteness factor of the typical density and size of in- 
stantons in the vacuum, k = ni nst p 4 ~ 10~ 2 (see e.g. || 
for details), the inelastic cross section reads 

Om(s,i)~ C7rp 2 lns J <%j_ dq 2 ± K(gu_, Q2±,t) 

dQ 2 K 2 B{Q) (8) 

(qi±+q2±) 2 
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where K is a pertinent instanton form factor at fixed 
—t <C s, containing through-going partons in the form 
of Wilson lines, k appears squared since the cross sec- 
tion involves the squared amplitude. B(Q) is the partial 
multi-gluon cross section for fixed Q 2 -C s and C an 
overall constant which accounts for both the instanton 
and antiinstanton contributions to the initial state. To 
exponential accuracy 

« 2 B(Q)«Im J dTe QT - s ^ (9) 

where S(T) is the effective action describing instanton- 
antiinstanton interaction for a time-separation T, which 
is defined to reduce to twice the free instanton action at 
large T. This effective action, also known as the "holy 
grail function" , is supposed to sum up contributions of 
any number of produced gluons. For small Q or large T, 
the dipole approximation is valid and (^) rises exponen- 
tially with Q B. However, as emphasized first by Za- 
kharov |^] the unitarity constraints on the partial cross 
section requires B(Q) to fall at large Q. Shifman and 
Maggiore argued that the unitarization could be quali- 
tatively enforced by resumming chains of instantons and 
antiinstantons. In Js) we followed this idea and indeed 
found that the dominant contribution to (^) occurs at 
the sphaleron point for which B(Afs) 1/k. Hence, 

c-in(M) « Cnp 2 Kins J dqu_dq 2 j_ K(qu_,q 2 A_,t) . (10) 

The rise in the partial inelastic cross section due to the 
production of the Yang-Mills (YM) sphaleron PI results 
into an increase of the inelastic cross section by one power 
of the diluteness factor k, or about a 100-fold increase ||. 

This qualitative solution of the problem implies how- 
ever that the maximal cross section does not corre- 
spond to a well-separated instanton-antiinstanton pair, 
but rather to a close pair with T sa p whereby half of the 
asymptotic action is annihilated. Detailed studies of close 
instanton-antiinstanton configurations of such kind have 
been made recently in |ic| l . An important new element of 
that paper was the emphasis on the t=0 3-plane, treated 
as a unitarity cut and explicitly describing the turning 
(escaping to Minkowski space-time) gauge field configu- 
ration. Although this paper was not dealing with the 
cross section calculation, following its logic one should 
e.g. modify the form-factor K in (@) by using Wilson 
lines in a half-instanton-antiinstanton field as described 
for instance by the Yung ansatz. 

However, even this solution of the problem would not 
be ideal, as it treats independently two small semiclassi- 
cal factors, the instanton-antiinstanton interaction S in 



1 Of course, there is also production of the YM antisphaleron, 
which carries opposite Chern-Simons number. 



the vacuum and the form- factor K. The idea behind the 
Landau method is to treat these semiclassical contribu- 
tions simultaneously. The essentials of the method are 
as follows (more details can be found in Q): In quantum 
mechanics the overlap between the ground state and a 
highly excited state can be rewritten as a difference of 
two shortened actions for a pair of paths with energy 
and Q respectively. The corresponding integrals inter- 
polate between the turning points and the singularity at 
infinity. The location of the singularity of the gauge con- 
figuration plays the role of such an infinity. Specific Eu- 
clidean paths, used originally in pl[,p| , have singularities 
in the field located at r=0 and time T-T/2, see figj^. 
Outside the region marked by the dashed lines the so- 
lution is the universal singular instanton, describing the 
ground state. Between the dashed lines it is supposed to 
be the (so far missing) energy-Q solution: both solutions 
join smoothly at the dashed lines. Our aim is to find the 
interpolating solution, at least approximately, and inves- 
tigate at t = the internal field structure of relevant 
turning states. Their subsequent evolution will eventu- 
ally lead to predictions of what is actually produced in 
the collision. 
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FIG. 1. Euclidean pace-time (r,t) plane. The circles in- 
dicate the location of the gauge singularities, and the verti- 
cal lines show where the 0-energy and Q-energy solution are 
joined. 

IV. MORE DETAILS ON SINGULAR 
YANG-MILLS 

A. Singular fields: \t\ > T/2 

The branch of the field that interpolates between the 
vacuum at t = — oo and the singularity at t = —T/2 
with zero energy and minimizes (^|) is an instanton. The 
conjugate branch is an an antiinstanton and interpolates 
between the singularity at t — +T/2 and the vacuum at 
t = +oo. In covariant gauge, both branches are given by 

A«(x(±)) = 2 v ^(±)x„(±Mx 2 (±)) 
_ 2r 1 l u {±)x v {±)p 2 

{p 2 ~x 2 {±))x 2 {±) [ ' 
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where the — field refers to t < —T/2 time-branch, and 
the + field refers to the t > T/2 time-branch, with 



x{±)„ = {x,t ± ) v = {x,tTT/2±p)„ 



(12) 



The t'Hooft symbols are identified as r/(+) = 77 (antiin- 
stanton) and rj(— ) = rj (instanton). For \t\ > T/2 ( |ll| ) 
are singular self-dual solutions to the QCD Yang-Mills 
equations with zero-energy. Note that the singularity in 
( |ll| ) stems from the change p — > ip in the self-dual 0(4) 
instanton. A(+) is the time-conjugate of A{— ). 

The axial gauge A4 = is commensurate with the 0(3) 
symmetry, and the results for the ± branches follow by 
using the hedgehog gauge transformation 



U(x, t) = exp 



dt'<5>{x 2 +t' ) 



(13) 



modulo static gauge transformations. Under the action 
of (13) the A± part in the Lorentz gauge ( pi] ) is gauged 
to zero. The residual static gauge transformations are 
fixed by fixing the positions of the singularities in the 
axial-gauge to coincide with those in the Lorentz gauge 
at t = ^fT/2. In particular, 



A?(x ) -T/2) = At(x,+T/2) 



2p 2 



x 2 (x 2 + p 2 ) 



(14) 



B. Singular fields: \t\ < T/2 

The gauge configuration in the time interval \t\ < T/2, 
follows from the YM equations using (|l4|) as the singular 
boundary conditions. They are no longer constrained by 
self-duality, and hence carry finite energy Q. Fixed Q 
relates to fixed T through Q = dS/dT where S is the 
Yang-Mills action in the time interval \t\ < T/2. 



1. Above the Sphaleron 

For small Euclidean times T/p < 1 or large energy, 
the singular boundary conditions ([l3|) common to both 
Lorentz and axial gauge, imply for the axial gauge de- 
composition (||) that in this limit B — > C A through- 
out [§. In this limit, we will call D = 2B = 2C the 
action (||) reduces to 



4a 



dt 



\t\<T/2 



dx 



D 



■S2 



8a; 2 



D 1 



the energy is 
M s ~ 



3 

4a 



dx 



D 2 



8x 2 



D z 



(15) 



(16) 



1 f°° D 3 
N w — / dx 



2- 



8x ' 



(17) 



with the boundary condition D(r, ±T/2) = —Ap/r. The 
action (1151) is extremal for 



D(r, -T/2) 
2r/(t + T/2) 



D(r,t)/D(r,-T/2) 



dx 



(18) 



The transcendental equation (|l^) can be solved numer- 
ically. The solution is shown in Fig. 1 (thick line) for 
p/T - 10. 

A good approximation at the escape point is 



D(r, 0) 



Apr 



" 2 + ^pT 



(19) 



which interpolates exactly between the asymptotics of 
the transcendental solution with K = 1.854. ( jTgj ) is also 
shown in Fig. 1 (thin line). Its corresponding initial ra- 
dial density is 



e 00 (r,o) 



4tt 



24p 4 r 



4 _2 



g 2 (r 2 + (M S /Q) 2/5 P 2 ) 4 



(20) 



which integrates to Q. Note that the tunneling duration 
T relates to the energy Q through 



T _ K (M s 



2/5 



(21) 



The Chern-Simons number at the turning point for (jli 
is 



2 \M S 



2/5 



(22) 



At the sphaleron point with Q = Ms the configura- 
tion ( |T9| ) carries Chern-Simons number N = 1/2: it is 
a sphaleron. 
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FIG. 2. D(r,0) for p/T = 10. See text. 



the Chern-Simons number is 
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For Q > Ms the initial configuration follows from the 
sphaleron by a simple rescaling of the size p and the en- 
ergy density, 



P 

a/5 



p/X 

- A 4 e 00 



(23) 



with A = (Q/M s ) L/b . In section VB we will evolve 
the gauge field configuration into Minkowski space using 
Luscher-Schechter (LS) solutions These solutions 

have a purely magnetic field configuration at t = and 
the (radial) energy profile 



9 00 (r,0) =47rr 2 e oO?,0) 



4tt 48ep 4 LS r 2 
7(r 2 + p 2 LS ) 4 



(24) 



Comparing with the rescalings ( [2 3D we see that we have 
to use the LS solutions with the parameters 



Pls 



A 
A 4 



(25) 



2. Below the Sphaleron 

Below the sphaleron the analytical analysis is more in- 
volved in general. For small energy Q or large times T, 
a perturbative expansion around the singular instanton- 
antiinstanton configuration has been carried out in ||. 
As we will argue below, the multiplicities below the 
sphaleron follow from the same rescaling ( p3|) with Q < 
M s . 



V. GLUONS IN/OUT 

In this section we estimate the number of incom- 
ing (virtual) and outgoing (real) gluons present in the 
semiclassical singular gauge configurations for arbitrary 
parton-parton center of mass energy Q. 



with, 



and 



Q^k) = —p 2 e ±k{T / 2 -^ 
9 



(27) 



4tt 2 p 2 
9 



,±k{T/2-p) 



(28) 



The transverse polarizations are denoted by A™ 1 . In terms 
of the Fourier components ( p8| ) the density of incoming 
transverse gluons is proportional to the occupation num- 
ber a a with 



(29) 



where 9 is a parameter fixed by requiring the total energy 
of the incoming gluons to match the energy Q Jll],|l2| . We 
note that T drops in the combination (|29|). Hence, the 
density of transverse gluons per unit wavenumber is 



n(fc) 



16tt 

a 



p(fcp) 3 e- 4fcp(Ms/2Q)1/S , 



(30) 



and the corresponding energy density per wavenumber is 



uj(k) = k n(k) 



16tt 



{kpfe 



4 -4kp(M s /2Q) 1/5 



(31) 



Under the rescaling (|23|) the energy density (|31j) of the 
incoming transverse gluons follows from the sphaleron 
point and integrates to Q. The virtual number of gluons 
stripped by a sphaleron is 



NUQ) - - 

mWJ 8a \M S 



4/5 



(32) 



and similarly for the antisphaleron. Each virtual gluon 
carries away about M5/3 « 2/3 GeV. We recall that 
these gluons are absorbed from the 2 eikonalized partons 
involved in the inelastic cross section (0). 



B. Outgoing Gluons 



A. Incoming Gluons 

The number of incoming gluons follow from the exact 
Euclidean solutions at large Euclidean times by expand- 
ing in powers of the space Fourier transform of the large 
Euclidean asymptotics of the singular fields ( |Tl| ) as in 
PIP^I . In Lorentz gauge, 



A%(t,k) = 2ik a Q(±,k) 
Af(t,k) = -2k a k t Q(±,k), 



A? 



fT(±,k)e 



-kt 



-hi 



(26) 



To assess the number of outgoing transverse gluons 
produced by the singular gauge configurations in the 
semiclassical approximation, we need to know the gauge 
configurations at the escape point and their further 
Minkowski time-evolution, much like the decay of the 
sphaleron in the standard model |14j]. The escaping 
sphaleron in Minkowski space is related to an analytical 
solution discovered by Luscher and Schechter [|l3| (LS) 
as discussed in jl0|,[l5|. What is remarkable in our case, 
is that through the scaling laws ( p3[ ) we have tied fea- 
tures of this solution (energy density and multiplicity) to 
those of the escaping singular Yang-Mills fields above the 
sphaleron point. 
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1. At the Sphaleron 

The LS solution in Minkowski space at the sphaleron 
point can be expressed in terms of elementary func- 
tions jl6j (some helpful relations can be found in the 
Appendices). It is strongly peaked around the light cone 
t ~ r as it travels luminally, and for large times t ~ r 3> p 
in covariant gauge it simplifies 



A%(t,v) 



EL 

9 
F' 



Af(t,v) ~ +—n a ni 

F' 

H — - (v (S ai - n a ni) + pe iaj rij) , 
9^ 

with v = r — t and 

2pf 

F (v) = -5 ~ 

/(u)= H 1- ^)) ' 



(33) 



(34) 



and £ = ta,n~ 1 (p/v). Note that /(±oo) = 1/2 - l/y/2 ~ 
-0.207, with /(0) > /(oo). Choosing F(0) = we see 
that F(v) is an odd function with asymptotics F(±co) — 
±0.216. We discuss some more details of the solution in 
Appendix A, and show that behavior of gauge invariants 
is the same as founding. 

The gluon number and spectrum is a subtle issue, and 
should be performed in "physical" gauges. In the covari- 
ant gauge, the large time asymptotics of the temporal 
and longitudinal gauge fields fulfill d ■ A = 0. They are 
proportional to F'(v) with support only on the light-cone. 
This is a gauge artifact in the covariant gauge, and can 
be removed by transferring to temporal gauge by using 
the hedgehog gauge transformation 



u(v) 



:F(v) 



(35) 



which yields Af t = 0. The temporal gauge is canonical 
in the sense that Gauss law is easily implemented by 
restricting to the transverse gluons, and the vacuum state 
is normalizable. In this gauge the large time asymptotics 
of the field is purely transverse and falls as 1/t ~ 1/r, 

A?(t, v) ~ -| (sin 2 F + y (pcos (2F) + vsin (2F))) e aij n 3 



1 

~~gi 



(-sin (2F) + F' (psin (2F) - vcos (2F)) {S ai - n a m) . 

(36) 



Note that the gauge transformation (35) modifies the 
Chern-Simons number ([?]). The transverse fields in both 
covariant and axial gauge weaken asymptotically as 1/t. 
Asymptotically the Yang-Mills solution originating from 
the sphaleron point Abelianizes, thereby allowing for a 



free wave interpretation. It is easier to carry the analysis 
in covariant gauge |. 

The large time transverse asymptotic of ( |33| ) fulfills 
trivially the covariant gauge condition, and admits a nor- 
mal mode decomposition in the form 



A?(t,k) 



(27T) 



2A- 



x (A™(fc)& am (fc)e^ fc * + A™(-fc)6 am *(-fc), 
with the A's as the 2 real polarizations, 
Xf(k)b am (k) 



+ikt 



(37) 



<7V7T k 

x(-ie a ijkj J(kp) + (5 a i - k a h) J'(kp)) , 



(38) 



and 



,7T/2 / 

J(fcp) = 2Re J dyeikpcotzny /j 



V2 



cosh(v / 2y) 



(39) 



The transcendental function ( J39| ) is not emmenable to 
closed form, but is well behaved for kp <C 1 



J(kp) ~ 7r — 4arctan ^tanh f ^J^ j 



J'(kp) ~ tt(V2- 1) 



(40) 
(41) 



In terms of the normal mode decomposition (|37|), the 
asymptotic density of transverse gluons is proportional to 
the occupation number |A • b\ 2 of the transverse modes, 



n(fc) =47rk 2 \X-b\ 2 



8p 2 k 



(j 2 (kp) + J' 2 (kp)) . (42) 



The number of gluons with small energy grows as kp, 
while the number of gluons with high energy falls as 
kpe~ 2kp . The total number of prompt gluons emitted 
by a sphaleron is 

Q POO -i i 

N out (M s ) = 4 / xdx ( j2 ( x ) + J ' 2 ( x )) = — ( 43 ) 
9 Jo a 

where the last integration has been performed numeri- 
cally^ The same number of gluons are produced through 



2 In temporal gauge there is a subtlety related to the constant 
modes which do not admit a spectral representation. Indeed, 
we have checked that the space Fourier component of j3^ ) 
exhibit a non-spectral term t5(k). 

3 The present spectrum is very similar to the one obtained 
numerically in |l(J using the temporal gauge. It is different 
from the one discussed analytically in [hoj using a different 
gauge where the gluon number was found to diverge logarith- 
mically at small k. The gluon number is also found to diverge 
in axial gauge in relation to the k — modes (see footnote 1). 
The results discussed above are free from gauge artifacts. 
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the antisphaleron. Each mode in the transverse asymp- 
totic ([57]) is normal, so that the energy density carried 
by these modes is co(k) = k n(fc), which integrates to give 
back the sphaleron mass by energy conservation, 



M s = 4~ f" dxx 2 (J 2 (x) + J'\x)) = — 
9 9 Jo a 9 



(44) 



The numerical result (^4j) is off by 7% in comparison to 
the exact sphaleron mass (Q) which is a measure of the 
onset of the asymptotic normal mode expansion (|37j). In 
Appendix C we show that the expanding sphaleron con- 
figuration is stable under pair production of light quarks 
and gluons. 



2. Away from the Sphaleron by Rescaling 

The escape configurations above the sphaleron follows 
from the gauge configuration (18-15^). The latter yields 
the energy density of the sphaleron after the rescaling 
( p3| ) and a Chern-Simons number of 1/2 at the sphaleron 
point. Since the electric field vanishes at the escape point, 
we conclude that it is very plausible that this gauge con- 
figuration is gauge equivalent to the LS gauge configura- 
tion with the parameters (p5|). 

The analysis of the preceding section may be per- 
formed just with the substitution of f(v) by 



l + V&.dn [ Jl + V^- 



for e < 1/2, and with 



1 + V2e ■ cn 



(8e)i£,m 



for e > 1/2, where 



1 



>2e 



2V2e 



(45) 



(46) 



(47) 



In Fig. 2 we show the multiplicity distributions for var- 
ious values of A. Numerically there is not much difference 
between the solutions obtained from the elliptic LS pro- 
files and an appropriate rescaling (^3|) of the sphaleron 
results. In particular, 



n(fc) ~ A 4 



./-' 



J 



(48) 



with A = (Q/Ms) 1 / 5 - The total number of prompt gluons 
emitted above the sphaleron is (in this approximation) 



iVout(Q) 

N out (M s ) 



Q_\ 

Ms) 



4/5 



(49) 



The ratio of in (virtual) to out (real) gluons per sphaleron 
is a number independent of Q: Ni n /N out ~ 2. Prompt 



inelastic scattering in QCD is from few-to-few (small Q) 
or many-to-many gluons (large Q). 




2 4 6 8 10 

FIG. 3. Density of emitted gluons per sphaleron (multi- 
plied by a/p) for A = 0.5, 1, 2 versus kp. 

Below the sphaleron barrier, there are two turning 
points: t = —Ty/2 with zero Chern-Simons number 
N = 0, and t = with unit Chern-Simons number 
N = 1. The two gauge configurations and thereby mul- 
tiplicities are related by the gauge transformation (|13|) . 
The remarkable similarity between the scaling law ( |49| ) 
for the outgoing gluons above the sphaleron point and 
the scaling law ( |32| ) for the incoming gluons both above 
and below the sphaleron point leads us to conjecture that 
holds for the outgoing gluons below the sphaleron 
point as well. 

The total number of prompt gluons emitted by the 
escaping singular Yang-Mills configurations below the 
sphaleron follows the scaling law ([49]), which is seen to 
vanish at the instanton point. This result follows from 
a saturation of the partial cross section via classical and 
singular solutions to the Yang-Mills equations, and is dif- 
ferent from the one derived recently using a minimization 
of the energy at the escape point by constraining the size 
and Chern-Simons number at the escape point ||10(] . The 
latter is likely to provide a lower bound on the partial 
cross section, while the former saturates it. 



C. Averaging Gluons 

In so far, we have considered the production of prompt 
gluons for fixed Q in the inelastic production cross sec- 
tion given by (j|). For the singular Yang-Mills solutions 
considered here, the partial cross section cr(Q) has been 
derived for small and large Q in || . In units of the scaled 
energy x = Q/M s , their result, with exponential accu- 
racy (and not too close to the x=l point) is 



a(x) = a + (x)9(x - l) + a-(x)6(l - x) 



with 



a±(x) 



(50) 



(51) 



7 



The (known part of the) 'holy-grail' function reads ^ 

F+(x) = -0.482 x 3/5 
F-(x) = -1 + 0.6185 x 4/5 

-0.0710 £ 6/5 + 0.0122 x 8 / 5 (-In x + const) . (52) 

The first contributions in F± are from the semiclassical 
singular gauge configurations alone, while the last two 
contributions in F_ are from the one-loop and two-loop 
contributions respectivel y [pi . The initial increase in the 
partial cross section in ( |5C| ) follows from the rapid in- 
crease in the tunneling rate at the expense of the de- 
crease in the matrix-element overlap. At the sphaleron 
point x = 1, F+(l) ~ F_(0)/2 which is about half the 
instanton suppression factor, 



o+(l) ~ \/o--(0) ~ \/e~^l a 



(53) 



in agreement with the unitarization arguments in |plg|. 
The final and rapid decrease in the cross section passed 
the sphaleron point is caused by the decrease in the over- 
lap between the initial and final states of the inelastic 
collision process. 

In this subsection we would like to show that in the 
inclusive gluon production a sharp peak remains, even 
after all gluon multiplicity factors are included. In Fig. 4 
we plot x 4 / 5 o~(x), with the typical instanton action in the 
QCD vacuum set to be 2n/a s (p) = 12. From this figure 
one can see that the r.h.s. of the peak is larger than 
the l.h.s.: it looks like Breit-Wigner peaks distorted by 
multibody phase space. 

Using (gOJ) we can assess the averaged number of 
prompt gluons produced in a parton-parton scattering 
process at large y/s, 

as a function of the invariant mass transfer in the t- 
channel, with the measure n(x, t) fixed by the partial 
cross section 



fi(x,t) = cr(x)/ / dx 2 a(x) 



(54) 



^Vin,out are the number of gluons at the sphaleron mass. 
The results for the ratio (pq ) are displayed in Fig. 5 versus 
—t/Mg. The undetermined constant in ( p2[ ) was adjusted 
to 2.418 to insure a smooth transition at x = 1. For \f—i 
within Ms, the averaged multiplicity is increased by only 
a factor about 1.1 compared to the multiplicity at the 
sphaleron mass. 



in, out 



(*) 



iVi, 



t/M§ 



dx 2 fi(x,t) x 4/5 



(55) 



4 The scaling laws derived above are specific to the energy 
density and the gluon multiplicities. They do not carry to 
the action den sity needed to assess semiclassically the partial 
cross section (pd). For that we need explicitly the escaping 



y 

6e-06 




0.6 0.8 



1.2 1.4 



FIG. 4. The inclusive gluon multiplicity x 4 ^ 5 a(x) versus 
energy in units of the sphaleron mass x = Q/M s . 



2.5r 



1.5 



0.5 



FIG. 5. The ratio (B5J) of prompt gluons emitted as a func- 



tion of -t/Ms- 



VI. CONCLUSIONS 



classical fields starting from (h 



for instance. 



We have presented a semiclassical assessment of the in- 
elastic cross section for multi-gluon production through 
instanton induced parton-parton scattering, using singu- 
lar classical solutions to the Euclidean Yang-Mills equa- 
tions. For energies Q much above the sphaleron mass we 
have found an approximate but explicit solution which al- 
lowed us to obtain the field configuration at the turning 
(escape) time t=0. The approximate solution turns out 
to be just a rescaled Yang-Mills sphaleron solution flio| |. 
We have solved the corresponding Yang-Mills equations 
describing its Minkowski explosion at later time, which 
is just a rescaled version of a Luscher-Schechter solution 
discussed in [Q. We have used the asymptotic form of 
the gauge field in physical gauge, to discuss the gluon 
number and spectrum. 
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Using an analytical continuation of the singular tails 
into Euclidean time, we have found that the same rescal- 
ing allows for the determination of the virtual in-gluon 
multiplicities both above and below the sphaleron (an- 
tisphaleron) point. Thus, the in-out multiplicities at 
the sphaleron and antisphalcron point, implies the in- 
out multiplicities away from this point for semiclassical 
parton-parton scattering. The scattering liberates about 
X.l/a s {p) prompt gluons by sphaleron (antisphaleron) 
produced. We recall that all the calculations are semi- 
classical, assuming a s (p) C 1 or large produced multi- 
plicities. The larger the longitudinal Q transferred in the 
CM frame, the smaller the size of the escaping sphaleron 
or antisphaleron and vice versa. The semiclassical pro- 
duction through singular QCD solutions in the semi-hard 
regime, may even extend to the hard regime through 
smaller size sphalerons. However, in this case the cross 
section is small. On the other hand, for typical instan- 
tons in the instanton vacuum relevant to the 'semi-hard' 
scale and possibly to the QCD pomeron problem, this 
number turns out to be only about 3 gluons per cluster 5 . 
In QCD there would be also quark production which we 
will discuss elsewhere. 

This mechanism of coherent multi-gluon production is 
additional to perturbative BFKL ladders and/or virtual 
gluon materialization in the color-glass approach [tt7l . 
The difference between this mechanism and others is seen 
e.g. in the fact that the released semiclassical gluons form 
thin shells of strong coherent field, which carry topo- 
logical features inherited from the topological tunneling. 
This difference is important for many applications, in- 
cluding the production of quarks as we discuss next. 



-eA? = e ian ^(l + 4> 2 ) 



r 

X a Xi 



X a Xi\ (pi 



(A2) 



The coefficient functions for the LS solutions continued 
to Euclidean space are 



A Q = -4g L 57 rt 



Ai = -4g L s7 : 



1 + r 2 - t 2 



1 - r 2 - t 2 



h = -^qLSl 2 r 2 - 1 



7 2 = 



(1 -r 2 -t 2 ) 2 + 4r 2 



(A3) 
(A4) 

(A5) 
(A6) 

(A7) 



and c[ls is a function of arctanh ^ 1+ ?l +t 'i j ■ 

In terms of these quantities the electric and magnetic 
fields are given by formula^ (7.27) in ji6fl . We perform the 
calculations for the sphaleron point which corresponds to 



V2 



cos ( v^arctanh 



21 



l+r 2 +t 2 



(A8) 



We rotate these formulas back to Minkowski space us- 
ing t — > it. The solutions are again concentrated on the 
light cone. We perform then the limit f -> oo keeping 
v = r — t fixed. The results for the 'electric' coefficients 
of formula (7.27) in Jl6| are 
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APPENDIX A: LS SOLUTION 

In this section we give a brief characterization of the 
LS gauge configuration. The reader is also referred to 
Ml for more details. For the Euclidean Witten ansatz 



eA° 



-An 



(Al) 



5 There are of course also instantons with smaller size, larger 
action and outgoing multiplicity. However in order to observe 
those the price to pay is exponential in action. Recall that 
electroweak sphaleron releases about 50 W,Z,H, which would 
be quite spectacular, but unobservable for the same reason. 



D cj)i/r 



v + 2sech 2 u(sinhw — v) 



t(i + v 2 y 







1 — 2sech 2 u(wsinhw + 1) 



t(l+v 2 ) 



2\2 



(A9) 
(A10) 
(All) 



where u = v / 2arctanh(l/w). The 'magnetic' coefficients 
are 



1 — 2sech u(v sinhw + 1) , , . 
i — (A12) 



t(l + v 2 ) 2 



i.o 

t 



. v + 2sech 2 «(sinh u — v) 
'' t(l + v 2 ) 2 



(A13) 
(A14) 



For large times t, the electric and magnetic fields are 
equal 



In the original formula (7.27) the coefficient of Foi is off by 
a factor of r 2 . 
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E 2 (t,v)~ B 2 (t,v)~ - 



1 



t 2 (i + v 2 f 



(A15) 



and sums up to Mg. The same result was obtained 
in@. 

The 'antisphaleron' follows from the sphaleron by sub- 
stituting 



A 

4>2 



-A 
-Ax 

-4>i 



(A16) 
(A17) 
(A18) 
(A19) 



in the Witten ansatz above, and yields a solution of the 
YM equation of motion. The same transformation makes 
a switch between the instanton and antiinstanton in the 
Witten ansatz. 



APPENDIX B: Im S FOR LS SOLUTION 

In this appendix we will show that the imaginary part 
of the classical action along the deformed contour |ll| 
does not depend on the energy of the LS solution. In 
other words, the instanton-antiinstanton suppression fac- 
tor persits all the way to the sphaleron point in the de- 
formed contour approach suggested in |lq| . Indeed, for 
e < 1/2 the solution is given by the function 



q = q_sn(q + <j) + K, k) 
where we use the same notation as in 111 



q ± = A/ 1 ± 



k = 



and 



arccoth 



1 + r 2 + t z 
27 



The suppression factor is given by 

Im S = 7t— [ 

9 Jo 

where 



dr Im Res 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



The residues are calculated at the poles of (Bl) inside the 



contour. One can show that £ is just equal to (q — 1) — e. 
We will now use the Laurent expansion of the elliptic sine 



sn(it) 



u 



+ (3(u - u Q ) 



(B6) 



Using the differential equation for sn(w) we obtain the 
relations a 2 k 2 = 1 and = f (1 + k 2 ). We thus have to 
compute 



Res (cos 4 w£ ) = a 4 fc 4 Rcs 



cos 4 w 



(4a 3 /3fc 2 <£ _ 2a 2 fc 2 )Res ( -^-^_) (B7 ) 



The coefficients of the residues are energy independent 
and equal to 1 and —2/3 respectively. Once we use 



cos 2 w 



At 1 



(1 +r 2 ) 2 + 2r 2 (r 2 - 1) + r 4 



(B8) 



and the r dependence of the poles r(r) = q + 
^/ q 2 — 1 — r 2 where q is an e dependent constant, we 
obtain 



dr ( cos 4 w \ i 1 + q 2 

1 



n r 
00 dr 



Res 



(i 



= i- 



l-q 2 



Finally we get 



Im S 



87T 2 

9 2 



(B9) 
(BIO) 

(BH) 



Note that the whole energy dependence through e has 
cancelled out. 



APPENDIX C: STABILITY UNDER 
PERTURB ATIVE PAIR PRODUCTION 

Are the escaping Yang-Mills fields stable under pair 
production, at least perturbatively? In general, they in- 
clude electric as well as magnetic fields: but the answer 
can still be given by assessing the gluon/quark pair emis- 
sion rate in a general time dependent background to first 
order in a. Let T{A) be the on-shell matrix element be- 
tween a particle and an antiparticle in the external gauge 
field. Typically, the number of pairs per unit 4-volume is 



Tr(a;| In (l - T(A) P(+) T(A) P(-)) \x) (CI) 



d A x 

where P(±) are the on-shell projectors on particles and 
antiparticles, and the trace is over color and spin. In the 
momentum representation 



P(±) = 2tt {p + m F ) 6(±P ) 8{P 2 



(C2) 



for a flavor of mass mp. Since the pair creation follows 
from large times (weak fields), we may approximate T{A) 
by its leading order 



T{A) 



(C3) 



Inserting (|C3j) into (|Clj), expanding the logarithm, inte- 
grating over all spacetime, and insisting on gauge invari- 
ance yield the pair rate 
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4a 



dNa, 



d i q 12 (2tt) 



|E-(q)| 2 -|B-( g )| 2 ) 



x% 2 - 4m 2 ,) (1 - Aml/q 2 ) 1 ' 2 ^ + 2m 2 F /q 2 ) . 

(C4) 

where (E , B) are the Fourier transform of the chromo- 
electric and -magnetic fields in Minkowski space 



(E, B) (q) = J d 4 x e iq - x (E, B) (a;) 



We have highlighted 4 in (C4) to show 



4=-x2x2x2 
2 



(C5) 



(C6) 



for g/\/2 fundamental quark charge, 2 flavors, 2 spins 
and 2 particle-antiparticle respectively. Notice that the 
production in (C4) is time like, for which q = (uj,0) is 
an allowed frame. Such frames support zero magnetic 
fields, indicating that the pair production mechanism is 
electric in nature. It is of or der a° in the strong coupling 
constant. In deriving (C5) we have ignored the back- 



reaction of the quarks on the YM fields. 

Since the uu and dd pairs are light, we may set rap ~ 
in flC4p . The total number of light quark pairs emitted 
by the classical field is 



AT. 



a 



uu+dd 



d\9{q 2 ) (|E?( ? )| 2 -|B?(</)| a ) 



(C7) 



The present arguments can also be used to derive a sim- 
ilar expression for the number of gluon pairs emitted by 
the classical field in the weak field limit. For N c = 3, the 
gluons can be organized in 3 conjugate pairs, such as: 
W + W~, K+K-, K°K° by analogy with the charged 
octet mesons. Say we choose the external field to be K , 
then K° may decay into the two charged modes W + W~ 
and K + K~ . Using the background field method, the 
result for the two gluon multiplicity is 



N a 



3a 



12(2tt) 4 



d\9{q 2 ) (|E«( 9 )| 2 -|B«( g )| 2 



(C8) 



A% 2 ) (|E-(g)| 2 -|B^( (Z )| 2 )=0 



(Cll) 



The expanding sphaleron starts magnetic and remains 
magnetic-like throughout. It is stable under light-pairs 
quantum emission. 
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x 2 x 2 



(C9) 



for 5\/3/2 charge, 2 deca y mo des and 2 physical helicitics 
respectively. Combining (C7) with (|C8|) we find that the 



light quark and gluon multiplicities are related 

2N F 



N- 



-N„ 



uu-\-dd g 

Using the results of appendix B, we find 



(C10) 
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